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Abstract 



We compute the vacuum polarisation on the lattice in quenched QCD using non- 
perturbatively improved Wilson fermions. Above of about 2 GeV^ the results 
are very close to the predictions of perturbative QCD. Below this scale we see signs 
of non-perturbative effects which we can describe by the use of dispersion relations. 
We use our results to estimate the light quark contribution to the muon's anomalous 
magnetic moment. We find the result 446(23) x 10^^'^, where the error only includes 
statistical uncertainties. Finally we make some comments on the applicability of the 
Operator Product Expansion to our data. 
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1 Introduction 



The vacuum polarisation Il{Q^) provides valuable information on the in- 
terface between perturbative and non-perturbative physics. It has been the 
subject of intensive discussions in the literature. 

The vacuum polarisation tensor is responsible for the running of ciem, which 
must be known very accurately for high-precision electro-magnetic calcula- 
tions. To calculate the hadronic contribution to the anomalous magnetic mo- 
ment of the muon we need to know the vacuum polarisation at scales from 
~ 100 MeV to many GeV. Perturbative QCD will be unreliable at the low end 
of this scale, so a non-perturbative calculation on the lattice would be useful. 

The vacuum polarisation n((5^) is defined by 

n^.(9) = ^Jd^x e"^-{0\TJ^{x)MO)\0) = {q^q, - q^g^,) Il{Q^), (1) 

where is the hadronic electromagnetic current 

- 2_ 1 - 

Jn{x) = 2^ e/ i)f{x)-t^'tpf{x) = -u{x)-tf,u{x) - -d{x)jf,d{x) + ■■■ (2) 

and = —q'^ (so that > for spacelike momenta, < for timelike), 
n can be computed on the lattice for spacelike momenta > 0. 

n can also be calculated in perturbation theory. 11 has to be additively 
renormalised, even the one-loop diagram (with no gluons involved) is loga- 
rithmically divergent. This renormalisation implies that the value of 11 can be 
shifted up and down by a constant depending on scheme and scale without any 
physical effects. However the dependence oiIl{Q^) is physically meaningful, 
and it must be independent of renormalisation scheme or regularisation. 

Experimentally 11 can be calculated from data for the total cross section 
of e+e~ hadrons with theoretical predictions of QCD by means of the 
dispersion relations [1] 

where 

7-)/ \ CTe+e- ^hadrons (s) I 3s \ , , 

The first derivative of the vacuum polarisation (the n — 1 term in eq.(3)) is 
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Fig. 1. A typical Feynman diagram contributing to the Cjj part of the vacuum 
polarisation n is shown as diagram (a). A diagram contributing to the An part of 
n is shown in (b). In this paper we only consider diagrams of the first type, with 
both photons attached to the same quark line. 

referred to as the Adler D-function [2] : 

D{Q^) = -Utt^Q^ dU{Q^) /dQ\ (5) 

The anomalous magnetic moment of the muon {g — 2)^ can be calculated 
to very high order in QED (5 loop), and measured very precisely, {g — 2)^ is 
more sensitive to high-energy physics than {g — 2)e, by a factor mjj^/ml, so it 
is a more promising place to look for signs of new physics, but to identify new 
physics we need to know the conventional contributions very accurately. QED 
perturbative calculations take good account of muon and electron loops, but at 
the two-loop level quarks can be produced, which in turn will produce gluons. 
The dominant contribution comes from photons with virtualities ~ m^, which 
is a region where QCD perturbation theory will not work well. 

n((5^) can be split into two contributions with a different dependency on 
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the quark charges, as shown in Fig.l, 

-127r^n(Q2) =J2e} C^(/x^ g^ mj) + ^ ejej, g^ m^, mj,). (6) 

/ /,/' 

The Cu term begins with a tree-level term which is 0{a^), while the first 
contribution to the An term is 0{ag). Furthermore, if flavour SU (3) is a good 
symmetry, the contribution to Au from the three light flavours {u, d, s) cancels 
because e„ + + = 0, and the only surviving contributions to Au come 
when both / and /' are heavy quarks (c, b, t). In this paper we will concentrate 
on the Cu term, both because it is larger, and because it is much easier to 
measure on the lattice. 

For large (spacelike) momenta Cn(/i^, Q^, "^/) can be expressed by means 
of the Operator Product Expansion, OPE 



Cuifi', g^ ruf) = co(/.^ Q^ m^) + ^^1^^ ruf {i^fi^f) (7) 

where Cq, cf , c'^^ and are the Wilson coefficients and n is the rcnormali- 
sation scale parameter in some renormalisation scheme such as MS. The /' 
sum extends over the flavours of the sea quarks (internal quark loops not di- 
rectly connected to the photon lines). The Wilson coefficients can be computed 
in perturbation theory, while the non-perturbative physics is encoded in the 
condensates m/ (■0/'0/(/i)), q;s(G^j, (//)), etc. 

Perturbatively the functions D, cq.(5), and R, cq.(4), are known to four 
loops for massless quarks [3,4], while the coefficient cq is known to three loops 
in the massive case [5]. The coefficients cf, c'^^ and c'^ are known up to 
0{al) [6,7,8] for massless quarks. 

In this paper we shall compute n((5^) on the lattice and compare the result 
with current phenomenology [9]. Preliminary results of this calculation were 
presented in [10]. 

The structure of this paper is as follows. After this Introduction we discuss 
the lattice setup in Sections 2 and Appendix A. The results are presented in 

Section 3 and Appendix C. In Section 4 and Appendix B we compare with 
perturbation theory. In Sections 5 and 6 we present a simple model which 
describes our lattice data well. In Section 7 we use this model to give a lat- 
tice estimate of the hadronic contribution to the muon's anomalous magnetic 
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moment. Finally in Section 8 we make some comments on the applicability of 
the Operator Product Expansion, and give our conclusions in Section 9. 



2 Lattice Setup 



We work with non-perturbatively improved Wilson fermions. For the action 
and computational details, as well as for results of hadron masses and decay 
constants, see [11,12,13]. Full details of the vacuum polarisation calculation 
can be found in [14]. 

To discretise the hadronic electromagnetic current 

jri^) = Y.^fJl(^) (8) 
/ 

we use the conserved vector current 

Jii^) =\ii^f{x + (^^){'^+lf^)Ul{x)i,f{x) 

where a is the lattice spacing. Prom now on we work with a Euclidean metric 
and write 

Ji{i)-Y.^''^^^'''''^Ji{^) (10) 

X 

with 

9. = ^in(^). (11) 



On the lattice the photon self-energy 11 is not simply given by {Jj^{x)Jl{0)). 
This is because the lattice Peynman rules include vertices where any number 
of photons couple to a quark line (not just the single-photon vertex of the 
continuum) [15]. We are therefore led to define 



n^.(g) = n«(g) + ng)(g), (12) 
with 



and 

^^^Kq) = -aS,AJ^^\0)) . (14) 
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where 

J^^\^) =l{^P{x + aft){l + ^^)Ul{x)^|;{x) 
+ 'iij{x){l - -ff,)Uf,{x)ilj{x + ajl)) 

It then foUows that (12) fulfils the Ward identity 

The conserved vector current (9) is on-shell as well as off-shell improved [16] 
in forward matrix elements. In non-forward matrix elements, such as the vac- 
uum polarisation, it needs to be further improved, 

j;""P(x) = J^{x) + ]^iaccvcd^.[i^{x)<T^,'il^{x)) , (17) 

where we have some freedom in choosing on the lattice. Any choice must 
preserve (16), and ideally it should keep 0{a?) corrections small. Considering 
the fact that the conserved vector current J^{x) 'hves' at x + ajl/2, the natural 
(or naive) choice for the improved operator would be 

+T^,^{x + afj, + au) - T^„{x - av) - T^„{x + ap, - ai>)^ 
= M^) + lccvccos{aq^/2) sin(ag^) ^ e*^%^(a;) , (18) 

where T^vlx) = ^(a;)c7^,y^(a;). Unfortunately it turns out that this choice 
introduces very large 0{a'^Q'^) errors into 11. We found that 

J7P(g) = J^{q) + ccvcsm{aq,/2)J2e''''T,4x) (19) 

X 

was a better choice of improved current, because it makes the 0{a^Q^) terms 
much smaller, and so we will adopt this definition. The difference between 
the definitions (18) and (19) is 0(a^(5^), so we are free to choose whichever 
definition leads to the largest scahng region in Q'^ (both agree at small Q^). 

In Fig. 2 we show the unimproved polarisation tensor along with the two 
different choices of improvement term in the case of free fermions. The best 
agreement with continuum physics comes from using (19), which is the pre- 
scription we will use in the rest of this paper. 
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Fig. 2. The effects of improving the current operator. We compare various curves, 
all calculated in the free-field case at the one-loop level with q oc (0, 0, 0, 1) and 
K = 0.123 (corresponding to am = 0.065). The dotted curves show the continuum 
one-loop result, shifted vertically to match with lattice results. The dashed line 
shows the one- loop lattice result without any improvement of the current. The 
dot-dashed curve shows the result of improving with (18) and the solid curve the 
result of improving with (19). The upper two curves differ by 0{a), the difference 
~ amlna^q^. The lower lattice curves are both 0{a) improved, but we see that the 
naive improvement (dot-dashed curve) has very large O(a^) discretisation errors. 
Improving with (19) produces a lattice result much closer to the continuum result. 



In Appendix A we give explicit expressions for Il^^j (q) and H^^^ (q) in terms 
of the link variables and the quark propagators. 
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3 Lattice Calculation 



To facilitate the extrapolation to the chiral and continuum limits, we have 
made simulations at three different values of (5 with three or more different k 
values at each /3. The parameters are hsted in Table 1. The lattice data for 
the vacuum polarisation for the individual momenta and (5 and k values are 
given in Appendix C. 



/? 


n 


V 


# Conf. 


6.0 


0.1333 


16^ 


97 


6.0 


0.1339 


16^ 


44 


6.0 


0.1342 


16^ 


44 


6.0 


0.1345 


16^ 


44 


6.0 


0.1345 


32^ 


16 


6.2 


0.1344 


244 


51 


6.2 


0.1349 


244 


38 


6.2 


0.1352 


244 


51 


6.4 


0.1346 


32^ 


50 


6.4 


0.1350 


32^ 


29 


6.4 


0.1352 


32^ 


50 



Table 1 

Parameters of the lattice simulation. The improvement coefficient in the fermionic 
action was taken to be csw = 1.769 for (3 = 6.0, csw = 1-614 for /? = 6.2 and 
csw = 1-526 for P = 6.4 [17]. 



First, let us discuss the value we use for Ccvc- From the fcrmion-loop con- 
tribution to the gluon propagator, computed in [18] to 0{m) in lowest order 
of lattice perturbation theory, we obtain 



c''o\f, am) ^3(ln{a'^f) - 3am{l- ccvc) H^'^q'^) - 3-25275141(5) 
+ 1.19541770(1) ccvc - 7-06903716(4) c^cvc 
+ am (6-46270704(30) - 5-29413266(6) ccvc 

+ 1-67389761(2)02^^^7)) + 0{a\m^) , 



(20) 



where am — 1/2k — 1/2kc- We use Cq' to refer to the lowest order, g^, per- 
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Fig. 3. Chiral extrapolation of Cn(g^, m) at f3 = 6.0. The points shown range from 
a^q^ 3.0 to ^ 6.5. All are calculated with ccvc = 1- 

turbative contribution to Cq. As said earlier, the dependence of 11 and Cq 
is physical. Therefore in an 0(a)-improved calculation there should be no 
0{a) terms which depend on g^. On the other hand a constant added to Cq 
has no physical effect, so there is no objection to constant terms of 0{a) in 
eq.(20). We see that there is an unphysical amln(a^g^) term in (20) unless 
Ccvc = 1 + ^(fl'^)- In the following we take ccvc — 1 and make the ansatz 

Cn{f,am) = C^(g^am = 0) + amM{f) . (21) 

In Fig. 3 we show the quark mass dependence of Cn(g^,am) for several mo- 
menta, which justifies assuming a linear am dependence. In Fig. 4 wc show 
the slope M{<f ) with and without improvement of the conserved vector cur- 
rent. The derivative M should tend to a constant when Q"^ ^ rri^. With no 
improvement term (open points) we see that M has a logarithmic dependence 
on g^, corresponding to an unphysical amln(a^g^) term in Cq. We see that for 
Cava = 1 the slope is, within error bars, independent of down to small 
momenta. This shows that the choice ccvc = 1 has eliminated or greatly 
reduced the unphysical logarithm in cq. We conclude that the tree-level value 



9 



60 



M(q ) 



40 - 



20 - 




10 



10' 



2a2 

a q 

Fig. 4. The slope M{q^) for the improved (sohd points, ccvc = 1) and unimproved 

(open points, ccvc = 0) vector current as a function of a'^cf' aX (5 = 6.0. The dashed 
line shows the height of the plateau in one-loop lattice perturbation theory, setting 
ccvc = 1 in eq.(20). 



Ccvc = 1 is a good choice 



The lattice sizes in Table 1 were chosen such that the physical volume 
is approximately equal for all three (3 values. As we have done simulations at 
(3 = 6.0, K = 0.1345 on two different lattice volumes we check for finite volume 
effects. We have not found any, see Fig. 5. This figure also shows that on the 
larger lattice we can measure the vacuum polarisation at much lower values 
of Q^, which is an important advantage. The quark boundary conditions are 
antiperiodic in all four directions, while gluon and photon fields have periodic 
boundary conditions in all directions. 



If Ccvc should ever be computed non-pcrturbatively, or in perturbation theory, 
it should be kept in mind that we have used eq. (19) as the definition of the im- 
provement term. 
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Fig. 5. A check for finite volume effects. At /3 = 6.0, k = 0.1345 we have measured 
the polarisation tensor both on a 16'' lattice (white points) and on a 32'^ lattice (black 
points). The agreement is excellent, and we conclude that finite volume effects are 
negligible. 

4 Compcirison with perturbation theory 



The first thing to do is to compare lattice results with continuum perturba- 
tion theory, which we do in Fig. 6. For the perturbative contribution c^'^^^iq^, m) 
we use the renormalisation-group improved result given in eqs. (67) and (69) 
of Appendix B. We use Kjjg — 243 MeV [18,19] and ji — l/a, and we identify 
with q^. The Tq parameter is used to fix the scale [20], with tq = 0.5 fm. 

Cn calculated on the lattice and Cn in the continuum can differ by an inte- 
gration constant which can depend on and a. In lowest order perturbation 
theory this constant is found by comparing Cq^ in MS", eq. (55), with Cq^ cal- 
culated in lattice perturbation theory, eq. (20). Setting ccvc equal to 1 and 
ix = 1/a we find that Acq = cjf* - c^^ = -22.379. In next to leading order 
this becomes 

Aco = -22.379 - 12/11 ln(a*^/a''^*) + 0(a^'^*). (22) 
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Fig. 6. The deviation of lattice data from continuum perturbation theory. The data 
are at /3 = 6.0, k = 0.1345 on a 16^ lattice (white points) and on a 32^ lattice (black 
points) . Below ~ 2 GcV^ there is a visible deviation from perturbation theory. 
The dotted line shows the effect expected from a gluon condensate as expected in 
the OPE. 



In view of this result, the value Acq ~ —28.5 seen in Fig. 6 is reasonable. 

Fig. 6 shows that the lattice results deviate from perturbation theory at 
large Q^. This deviation, which sets in at a^q^ ~ 5, is probably a sign of O(a^) 
lattice artefacts. 

More interesting are the deviations from perturbation theory at low Q^. 
These are especially large at low quark mass. The OPE, eq.(7), would suggest 
that the effects of a gluon condensate should show up at small Q^. In Fig. 6 we 
show the OPE prediction for a gluon condensate ^{GG) = 0.012 GeV^. The 
deviations which we see in the data do not look like the effects expected from 
a gluon condensate. The sign is the opposite of what the OPE predicts, and 
the deviation is probably not growing as quickly as 1/Q^- 
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5 A simple model of the vacuum polarisation at low Q"^ 



We have seen that perturbation theory, even when supplemented with higher 
twist terms from the operator product expansion, has difficulty in explaining 
the low region of the data. Can we understand this region in some other 
way? 

The cross section ratio -R(s) is given by the cut in the vacuum polarisa- 
tion, or in other words there are dispersion relations which give the vacuum 
polarisation if we know i?(s). 





K 


amps 


amy 




6.0 


0.1333 


0.4122(9) 


0.5503(20) 


0.2055(14) 


6.0 


0.1339 


0.3381(15) 


0.5017(40) 


0.2217(15) 


6.0 


0.1342 


0.3017(13) 


0.4904(40) 


0.2293(14) 


6.0 


0.1345 


0.2561(15) 


0.4701(90) 


0.2387(90) 


6.2 


0.1344 


0.3034(6) 


0.4015(17) 


0.2210(20) 


6.2 


0.1349 


0.2431(6) 


0.3663(27) 


0.2403(24) 


6.2 


0.1352 


0.2005(9) 


0.3431(60) 


0.2474(47) 


6.4 


0.1346 


0.2402(8) 


0.3107(16) 


0.2252(21) 


6.4 


0.1350 


0.1933(7) 


0.2800(20) 


0.2423(19) 


6.4 


0.1352 


0.1661(10) 


0.2613(40) 


0.2448(35) 



Table 2 

Lattice data on the vector meson used as input for the dispersion relation fits [11,13]. 
Numbers in italics have been interpolated from nearby n values. 

One way of modelling 11 is to make a model for i?(s), and then calculate 
n from i?, using the dispersion relation eq.(3). This should be quite robust, 
since 11 gets contributions from a large range in s, little inaccuracies in the 
model R get washed out, and we can hope that even a crude representation 
of R will give a good result for 11. We will keep the model simple so that we 
can do all the integrals analytically. 

Perturbatively 

i?(s)=^e57Ve(l + ^ + ---) (23) 
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where A^c is the number of colours (3 in our case). Wc know that really the 
low s behaviour of R is more complicated than that, it is dominated by the 
p(770), a;(782) and 0(1020) mesons. Following [7], let us make the following 
model for R. We ignore the sphtting between the p and a;, which comes from 
the 74n-type diagrams which we have dropped. We also treat these mesons as 
narrow resonances, each contributing a S function to R. The continuum part 
of R takes a while to climb up to the value in eq.(23). We will represent this 
rise by a step function at some value sq- So, our model is 

^(«) = E 4 {A6{s - ml) + Se(s - So)) (24) 
/ 

where we would expect B to be slightly above 3 in order to match eq.(23). 



Using the dispersion relation, this i?(s) translates into a vacuum polarisation 

Cn(g') = B ln(a2Q2 + a'so) - >1/(Q' + m^) + X (25) 
where X is a constant which is not determined from the dispersion relation, 
and which never appears in any physical quantity. Again we identify the con- 
tinuum quantity with the lattice quantity g^. 



The constant A can be expressed in terms of the decay constant /y, which 
has been measured on the lattice. The cross section for the production of a 
narrow vector resonance, V ^ is [7] 

vriy 

The partial width Fy^e+e- is related to a meson decay constant gy [21,7] by 

( 47ra|^\ my 

r.^e.e- = \—^) , (27) 

where 

(0|J^-|T/,£) = £,^. (28) 

Here is the polarisation vector of the meson. On the lattice it is more natural 
to define decay constants jy in terms of currents with definite isospin [11] 

(0|V,|y,£) = £,^ (29) 



with 



for the (P and 



(30) 



V;=° = ^(iZ7/.t^ + d7/.^^) (31) 
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for the uj (ignoring any ss admixture in the uj). The relationship between the 
two definitions is 



l_ ^ e„ - ed J_ ^ ^ ^ /ooN 
J_ ^ eu + eg J_ ^ J_ 

In terms of these decay constants 



2 

R[s) = J2 127r25(s - m^) ^ + continuum (34) 
y 9v 

10 2c-/ 2\''^'p i^u — ^d)^ 

— IzTT d(s — ^ 



p 



+127r^5(s — rn^S}~^ ^"^ o"'^' — '~ continuum. (35) 

Neglecting annihilation diagrams implies that rriuj — rUp and /a, = fp. Ex- 
perimentally both relations are fairly accurate. The mass ratio m^/nip is 
1.02. = fp implies V^^_^f.+e- = |rp"^e+e-) while the experimental ratio 
is 0.089(5) [22]. Equating the mass and decay constant of the u! and p in 
eq. (35) gives 

2 

R{s) = 127r^(5(s - ml)^ (el + efj + continuum ; (36) 



so 

A = 1 9ir2 

fv 



12n^'-^. (37) 



6 Dispersion relation fits 



We first try making fits to the lattice data for the vacuum polarisation Cu 
using eq. (25) with B, K and Sq as free parameters. To avoid problems from 
lattice artefacts of 0{a?'(f) we have only used data with a^q^ < 5. We will call 
this simple ansatz Fit I. A, the weight of the vector meson contribution, is 
determined by eq. (37). The vector meson masses and decay constants which 
we use are shown in Table 2. They have been taken from [11,13]. We can 
compare the values for B and K with the one-loop lattice perturbation theory 
result, eq. (20), which gives B — and K — —27.38 -|- am 8.53. 
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/? 




V 


a^so 


B 


K 




6.0 


0.1333 


16^ 


0.395(32) 


3.13(6) 


-32.98(11) 


4.3 


6.0 


0.1339 


16^ 


0.397(43) 


3.13(8) 


-33.14(14) 


2.2 


6.0 


0.1342 


16^ 


0.391(39) 


3.10(8) 


-33.18(13) 


3.0 


6.0 


0.1345 


16^ 


0.409(61) 


3.11(9) 


-33.30(16) 


2.5 


6.0 


0.1345 


32^ 


0.403(79) 


3.12(10) 


-33.30(18) 


0.7 


6.2 


0.1344 


24^ 


0.282(21) 


3.10(5) 


-32.94(7) 


1.9 


6.2 


0.1349 


24^ 


0.278(24) 


3.07(5) 


-33.06(8) 


1.9 


6.2 


0.1352 


24^ 


0.253(25) 


3.06(5) 


-33.10(7) 


2.2 


6.4 


0.1346 


32^ 


0.164(13) 


3.06(3) 


-32.68(5) 


0.8 


6.4 


0.1350 


32^ 


0.150(14) 


3.03(4) 


-32.76(6) 


0.6 


6.4 


0.1352 


32^ 


0.124(12) 


3.02(3) 


-32.78(5) 


1.2 



Table 3 

The result of Fit I, eq. (25), including only data with a^q^ < 5. 

A typical result is shown in the top panel of Fig. 7. As can be seen from 
Table 3, the values for the fits are low in every case. Nevertheless we see 
that the data deviates quite strongly from the fit when a^g^ is large. The 
small deviations from the fit at smaller are not random - they occur in the 
same place in every data set, and depend on the direction of q. Points where 
q is near the diagonal direction (1, 1, 1, 1) lie lower than points measured for 
momenta away from this diagonal direction. 



To correct for this behaviour we add terms to our fit to parameterise 0(0^) 
lattice errors. There are two possible terms at O(a^): 

and a'^^ . (38) 

The second term has only cubic symmetry, and can give rise to a dependence 
of Cu on the direction of q, which would not occur in the continuum. To 
include these terms we fit with the ansatz (Fit II) 
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Fig. 7. A comparison between Fit I and Fit II. The data are at (3 = 6.0, n = 0.1333 
on a 16*^ lattice. Fit I only uses data with a^q^ < 5 (vertical dashed line). The 
second fit describes the data better. 
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Cuif) = B In(a2g2 ^ ^2^^^ _ _A_ + K + U^a^f + U2h{q) (39) 



where 



sin^ aq^ - \ (e^, sin^ aq^' 
a?q^ 



g2 4^ 



(40) 
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Fig. 8. The deviation of lattice data from the fit of Table 4. The data are at /3 = 6.0, 
K = 0.1345 on a 16^ lattice (white points) and on a 32^ lattice (black points). The 

fit describes the data extremely well. 

The function h{q) has been chosen so that it vanishes when goc(l,l,l,l), 
as most of our momenta are near this direction. Looking at the lower panel of 
Fig. 7 we see that the Ui term does a good job of describing the high q^ data, 
while the U2 term successfully describes the direction dependence of the data. 
Fits with the ansatz (39) show practically no deviation from the data. This 
can also be seen in Fig. 8, where we show the deviation in the case where we 
have data on two lattice sizes. 



In Fig. 9 we show the results converted into physical units. The agreement 
between the different (5 values is fair, and the value of Sq in agreement with 
phenomenology (for example, [23] finds Sq — 1-66(22) GeV^ in the 1=1 chan- 
nel). 

In Fig. 10 we show the values for the fit parameter B. The value we find is 
always very close to the tree-level value 3. 
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Table 4 

The result of Fit II, including all data. 
7 The anomalous magnetic moment of the muon 



With a good description of the vacuum polarisation tensor in the low 
region we can make statements about phenomenologically interesting quanti- 
ties such as the contribution of the d and s-quarks to the muon anomalous 
magnetic moment [24]. 

Traditionally the hadronic contribution to the muon's anomalous magnetic 
moment is found from R{s) via a dispersion relation (see [25] for a review) 

had = ^ r R^s) (41) 

^ 37r2 JAml s \mll ^ ' ^ ' 



where 




Io + (1 — x) s / mji^ 
The muon mass, m^, is 105.7 MeV. For us it is more useful to deform the 
contour and find this same number from the vacuum polarisation at spacelike 
momenta [24] 
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Fig. 9. The threshold sq (white points) compared with (black points). Squares 
are for /3 = 6.0, triangles for 6.2, and circles for 6.4. 



had ^em 



at. 



37r2 



where the kernel F is 



n(o) - n(g2) 

Cn(g',m/)-Cn(0,m/) 
+ annihilation , 



16 mf. 



mf 



(43) 



(44) 
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Fig. 10. The parameter B from Table 4. The /3 values are shown by the same 
symbols as in Fig.9. 

To calculate the value a^"'' we evaluate the integral 

I^i^f) = r i^] \Cu{Q', mf) - Cn(0, mA (45) 



for each of our data sets. We then need to extrapolate (or interpolate) to 
the chiral limit (for the u and d quarks) and to the strange quark mass. As 
can be seen in Fig. 11 the integral is dominated by Q"^ ~ 3m^, so we need to 
extrapolate in Q^, as our lowest measured value is at = 0.17 GeV^. We 
do our extrapolation of Cn by using the results of Fit II from the previous 
section. 

The results arc shown in Fig. 12, plotted against rn?pg (the square of the 
pseudoscalar mass). We see that / depends strongly on the quark mass, with 
heavier quarks making a smaller contribution (as one would expect). There is 
also some dependence on (3. To extrapolate to the continuum we fit the data 
with an ansatz of the form 

1^ = (^1 + A2a^) + (Si + B2a'')ml (46) 

This describes the data well (x^/dof = .53), and gives the continuum limit 
shown by the dashed line in Fig. 12. The physically relevant values of P are 
at the physical pion mass (M^ = 0.019 GeV^) for the u and d quarks, and at 
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Fig. 11. The product F{Q'^/ml) [C^(Q^m/) - Cn(0,m/)] for the case /3 = 6.0, 
K = 0.1345. The shaded region shows the uncertainty. 

the mass of a hypothetical ss pseudoscalar meson with = 2m\ — = 
0.468 GeV^ for P. The values at these points are 

r = 1'^ = 0.0389(21) , r = 0.0287(9) . (47) 

This gives as our final lattice estimate of a^'^'^ 

^had ^ ^ f^ju ^ Ijd ^ ljs\ ^ 445(23) X 10"^°. (48) 
Stt \ 9 9 9 / 

This error reflects the statistical errors of the lattice calculation and the ex- 
trapolations to the physical points. It does not include any estimate of the 
error due to the quenched approximation used in the calculation. This value is 
somewhat lower than the experimental value 692.4±5.9ea:p±2.4radiatTOe x 10"^°, 
found by applying eq. (41) to experimental R measurements [26,27]. Our value 
is similar to another lattice measurement [24], which finds 460(78) x 10~^°. 

The shortfall in the value of a^°'^ can probably be attributed to quenching. 
In particular a quenched calculation omits the process e~^e~ — >■ vrvr, which is 
the only contribution allowed at very low s. 
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Fig. 12. The results of the integral (45). The /3 values are shown by the same symbols 
as in Fig.9. The dashed line shows our extrapolation to the continuum limit, and 
the points shown with stars are our results for the light {u, d) and strange quarks. 

8 The applicability of the Operator Product Expansion 



What does the success of our fit function eq. (24) tell us about the applica- 
bility and usefulness of the OPE? 

To answer this question let us look at the large behaviour of the for- 
mula (25): 



Cn(Q') = 51n(g2) + X + 



Bsq-A Aml-^sl 

g2 + (Q2)2 

B 



( hi ^ 

B\n{Q') + K+ ^(-1)- [Aimlr-' - -s^ (g^)-^ 



(49) 



The expansion has the same form as the OPE, at least at leading order when 
there are no logarithmic corrections to the Wilson coefficients. We can now 
look at the higher twist terms in the expansion, and estimate how important 
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they are in comparison with the gluon condensate contribution. 



By looking at the first few terms we can relate our parameters Sq and A to 
the condensates in the OPE. In the chiral limit there is no l/Q^ term, so 

A = Bsq. (50) 

This says that the area under the meson (5-function is the same as the grey 
area in Fig. 13, a typical sum-rule style result. Substituting this into the l/Q^ 
term we get 

Bso(my^so^^-7r'^{GG) . (51) 

Again, a typical sum-rule result — if there were no gluon condensate, the 
meson would lie at |so, exactly in the middle of the gap. Putting = 
0.6 GeV^ (the physical value) and ^{GG) = 0.012 GeV^ gives sq = 1.26 GeV^ 



R(s) 



Sq s 

Fig. 13. Our simple model for the contribution of the u and d quarks to R{s). 

Next, in Fig. 14, wc plot a comparison of the threshold model, eq.(25), per- 
turbation theory (just 3lna^Q^ at this level) and perturbation theory plus 
the gluon condensate contribution, which scales like l/((5^)^. The physically 
irrelevant constant K has been set to in all cases. 



The curve from the threshold model looks physically sensible, it goes to a 
finite value at = 0, which is what must happen if R{s) doesn't extend all 
the way down to s = 0. 



However, even though we chose our parameters A and Sq so that the thresh- 
old model would match the gluon condensate prediction at high Q^, the two 
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10 ^ 10 ' 1 lo' 

Q^[GeV^] 

Fig. 14. A comparison of the threshold model, eq.(25) (solid line), perturbation 

theory (dotted line) and perturbation theory plus the gluon condensate contribution 
(dashed line). In this plot we have chosen the physically irrelevant constant K to 
be when all quantities are expressed in GeV. 



curves don't resemble each other closely. 

Let us now subtract out the perturbative piece, to see things more clearly. 
Fig. 15. We can only find a region where the gluon condensate region is impor- 
tant when we concentrate on the large region. If Fig. 15 is what the real 
world looks like, the prospects for getting at the gluon condensate look poor. 
The gluon condensate term is overwhelmed by higher order terms in the OPE 
before it has a chance to get significant. 

What is the conclusion of this exercise? In the region where we can easily see 
deviations from perturbation theory the OPE is not very useful, because many 
operators of very high dimension are all contributing, not just the leading 
1/{Q'^Y contribution coming from the gluon condensate. To see the gluon 
condensate uncontaminated by higher order operators we would need to look 
in the region ~ 10 GeV^ with very accurate data (error bars at least 2 
orders of magnitude smaller than in this paper). This is unfortunately not a 
realistic prospect. 
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Fig. 15. The same as Fig.l4, but with the perturbative piece subtracted. Note 
that in this figure we are concentrating on the behaviour at rather large where 
non-perturbative effects are smah. 



9 Conclusions 



We have computed the vacuum polarisation in the limit of two light flavours 
in quenched QCD for three /3 values (6.0, 6.2 and 6.4) and on lattices as large 
as 32^. It was important to improve the action and the vector current. We 
found good agreement with three-loop perturbation theory in the interval 
2 < < 20 GeVl The lattice data show some indication of non-perturbative 
effects at the lower end of the range. We can describe these very well with 
a model of Ris) which includes vector mesons and threshold effects. In order 
to make firm predictions, we need more data at small and at smaller lattice 
spacings with high statistics. 

From the low region of the vacuum polarisation we can extract a lattice 
value for c^^'^ ^ the hadronic contribution to the muon's anomalous magnetic 
moment. We find the value 446(23) x 10~^° which is of the right order of mag- 
nitude, though lower than the physical value. Our estimate could be improved 
by using a larger lattice size, enabling us to reach lower which would reduce 
uncertainties from extrapolation. Naturally, dynamical calculations would be 
very interesting. 
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Appendix A 



We denote the quark propagator from lattice points a; to y by G{x,y). For 
then obtain 




X 



,iq{x+ap,/2—av /2) 



Tr((l + 7,)t/t(0)75G'^(a: + a/i, 0)75(1 + ^^)Ul{x)G{x, av) 

- (1 - -i,)\JMl^G\^ + a/i, az>)75(l + ^;)Ul{x)G{x, 0) 

- (1 + 7,)[/J(0)75G't(x, 0)75(1 - ^;)U^{x)G{x + all, av) 
+ (1 - 7,)C/,(0)75G't(x, aj>)75(l - -i^,)lJ ^{x)G{x + a/i, 0)) 



— jccvc Qx 



X 




(52) 



- (1 - -f^)Ui,{x)G{x + afi, 0)a,x-f5G\x, 0)75) 

+ -jCCVC Qa 



X 



IV((1 + 7,)C/t(0)75G't(x, 0)75(7^.G'(x, au) 
- (1 - 7,)C/,(0)75G'^(x, ai>)^^a^G{x, 0)) 

^6 

^ X 

Tr{a^G(x, 0)a,xl5G\x, 0)75) , 
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where the improvement term has been integrated by parts. For Il^^^{q) we 
obtain 



ng](?) = ls,^Tr{{l + jMmiO, au) + (1 - 7.)C/.(0)G't(0, au)) . (53) 



To compute II^^J (g) and n|fj (g) wc have to do a minimum of five inversions for 
each gauge field configuration (and each k, value), which makes the calculation 
computationally quite expensive, in particular on 32'' lattices. 



Appendix B 



Perturbative results 



Before we describe the lattice calculation in detail, we present here the 
perturbative Wilson coefficients cq, cf and cf. We will work in the quenched 
approximation, in which contributions from sea quarks are neglected, and 
which corresponds to — 0. 



with Cp = 4/3 and Ca = 3 for SU{3). These coefficients can be found in 
eqs. (27)-(30) of [5] (recall that = -q^). In the WS scheme the coefficients 



We write 



TT 



(54) 



Cq , Cq , Cq aud Cq rcad 



4°^(/^^g^m) 




4^)(/.^g^m) 



9 



55 



4C(3)-ln 



g2 




4 



12 




(56) 



28 



c^,''(A.^<^^m) = -7 



_l|_!!c(3) + 10C(5) + il.^ 



72 
4to^ 



6 

1667 
~96~ 



8 

^C(3)-fc(5) 



(57) 



44215 
2592 

11 



22Vr'^^ 5 ^ 41 
^C(3)--C(5)--ln-^ 



11 



+ 24^^1^ + y^^^)^V 

4m2/1447 4^,^, 85^,^, 



185, g2 11 

In — T + 



24 



8 



In^ 



(58) 



where C(3) = 1.20206- • • and C(5) = 1.03693- • •, and rh refers to the quark 
mass at the scale ji in the MS scheme. In (55), (56) terms of 0{{m? /Q'^Y) 
modulo logarithms have been neglected, while in (57), (58) terms 0{{in^ /Q'^Y) 
are dropped. 

The Wilson coefficients multiplying the quark and gluon condensate are [6,8] 



cf(;.,,) = -12.^(2 + ^^ + ...), 

^ ^ ^ ^ (59) 



Note that both m{qq) and (cks/tt) (C^J) are renormalisation group invariants, 
which means that they do not depend on fi and the renormalisation scheme. 
The vacuum polarisation ll{Q^) itself is not an observable, but its derivatives 
are. Therefore the result (7) can only depend on /j, and the scheme in terms 
of an integration constant (independent of Q^). 
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Renormalisation Group Improvement 



We can use renormalisation group improvement to re-sum the logarithms in 
higher-order terms. This should lead to a significant improvement, since the 
fact that we are interested in measurements over a large range means that 
these logarithms are large. 



If wc calculate the Adler function in the chiral limit from eqs (55)- (58) we 
find the result 



3<1 + 



TT 



+ 



TT 



(60) 



We should be able to do better than this because the massless Adler function 
is known to four loops [3]. We can use this perturbative result for the Adler 
function to improve the result for cq in the chiral limit. In quenched SU{3) we 
have 



TT 

2\\3 



+ 



TT 



365 
'2A 



-11C(3) 



+ 



2\\ 3 



87029 1103^,^, 275^,^, 121 ^ 



'o:s{Q 



n 



(61) 



and 



365 



+(Ee/) 

/ 




87029 1103^,^, 275^,^, 



2 / Ois{Q 



(62) 



Note that although the first terms of it! and D coincide, the a\ term is differ- 
ent. The extra "vr^" term in R arises from analytic continuation, ln^(Q^/ yU^) 
[ln(s//i^) ± ivr]^, see [3]. The first part of (62), proportional to J2f^}, is the 
derivative of the Cn term in (6), while the second term, proportional to 
{Y^fCfY, comes from the derivative of the Au term. 
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Prom (62) we have a differential equation for cq: 



Q'^co(A^^ Q2 m= 0) = 3<' 1 + 



dQ 



+ 



TT 



365 



+ 



TT 



5Z2^_HEc(3) + HI5c(5) 

288 4 ^ 6 ^ ^. 



24 



11C(3) 
(63) 



We can solve this by using the known /3-f unction [28] 



Q 



2 d ( asiQ"^' 



ag2 1 TT 



'0Ls{Q 



2\\ ^+2 



with 



- 11 - 51 n _ 2857 

in the case of quenched SU{?>). 



TT 



(64) 



149753 891 , , , , 



The solution is 



Co(/x^Q^O) = 31ng2 _ i^ln(a,(Q2)) 



3403 



242 



+ 12C(3) 



TT 



(66) 



+ 



2301587 273,,^, .....x 
^^ + -C(3)-25C(6) 



TT 



+ C(q;^) + const. 



The constant of integration can be fixed by comparing with (54)-(58), giving 
the final renormalisation group improved result ^ 



Co(Ai',Q',0) = 31n%-^ln 



+ 



^j? 11 ya.s{tJ?) 

2301587 273 
15972 ^ ~Y 



+ 



3403 



242 

C(3) - 25C(5) 



+ 12C(3) 



asm 



TT 



'«s(g')\' 



TT 



(67) 



151 a,(//2) 
"^ + 484 ^^ + 
+0{al{Q%al{^,-)) 



566749 5 
383328 ^ 3^*^ 



Q;s(/^ ) 



2^\2 



TT 



where is the MS scale. 



^ Note the interesting result that due to a cancellation there is no ^(3) term in the 
coefficient of asip?) and no ^(5) in the coefficient of as(//^)^. This surprised us. 



31 



One can check that the differences between (67) and (54) are indeed 0{al) 
by making the substitution 



- Po 



TT 



in (67). 



TT 



TT 



In — - /?i 
//2 



2^\3 



TT 



+0i 



TT 



+ 0{a 



1 



(68) 



Renormalisation group improvement of the mass terms in (55)-(58) is sim- 
pler because there is no constant of integration involved. The result is 



co(Af^g^m) =co(Af^g^O) + 



+ 



19691 124,, , 1045,, , 



18 + 48 



TT 



+ ^_9-181n 



TT 



2 + 48C(3) + 66 In + 36 In^ ^ 



TT 



(69) 



What have we gained by using the renormalisation group? Most impor- 
tantly, the derivative of (67) is correct to four loops, one order better than 
(54). The uncertainties in (67) are 0{al{Q'^), q;^(/x^)) while in the original for- 
mula (54) terms of order 0(q;^(//^) ln^((5^//i^)) have been neglected. So, if we 
are interested in describing a large range of Q^, eq.(67) ought to be the better 
formula. In addition, (67) exhibits correct physics, because it is the sum of a 
piece depending only on and a piece depending only on /i'^, which is only 
approximately the case for (54). 



Appendix C 



In the following Tables we give the results for Cuiq, dfn) at the various /3 and 
K values. We label the lattice momenta by the vector n, = (2/a) sin(7rn^/L). 
The momenta were chosen close to the diagonal of the Brillouin zone to avoid 
large 0{a^) effects. 
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n 


K = 0.1333 


K = 0.1339 


K = 0.1342 


K = 0.1345 


(1,0,0,0) 


-38.15(11) 


-38.60(15) 


-38.84(14) 


-39.03(16) 


(1,1,0,0) 


-36.59(10) 


-36.89(13) 


-37.05(12) 


-37.17(13) 


(1,1,1,0) 


-35.52(09) 


-35.75(12) 


-35.89(11) 


-35.99(12) 


(1,1,1,1) 


-34.72(09) 


-34.91(12) 


-35.04(11) 


-35.12(12) 


(2,1,1,1) 


-32.89(08) 


-33.06(10) 


-33.15(09) 


-33.24(10) 


(2,2,1,1) 


-31.82(07) 


-31.98(10) 


-32.07(09) 


-32.15(10) 


(2,2,2,1) 


-31.07(07) 


-31.22(09) 


-31.30(08) 


-31.40(09) 


(2,2,2,2) 


-30.47(07) 


-30.62(09) 


-30.71(08) 


-30.81(09) 


(3,2,2,2) 


-29.61(06) 


-29.76(09) 


-29.85(08) 


-29.95(09) 


(3,3,2,2) 


-28.98(06) 


-29.13(08) 


-29.22(08) 


-29.32(08) 


(3,3,3,2) 


-28.47(06) 


-28.63(08) 


-28.72(07) 


-28.82(08) 


( 3 j3 j3 j3 ) 


-28.05(06) 


-28.20(08) 


-28.29(07) 


-28.39(08) 


(4,3,3,3) 


-27.54(05) 


-27.70(08) 


-27.79(07) 


-27.89(08) 


(4,4,3,3) 


-27.12(05) 


-27.28(07) 


-27.37(07) 


-27.47(08) 


(4,4,4,3) 


-26.76(05) 


-26.91(07) 


-27.00(07) 


-27.11(07) 


(4,4,4,4) 


-26.42(05) 


-26.58(07) 


-26.67(07) 


-26.78(07) 


(5,4,4,4) 


-26.11(05) 


-26.26(07) 


-26.36(06) 


-26.46(07) 


(5,5,4,4) 


-25.82(05) 


-25.97(07) 


-26.07(06) 


-26.17(07) 


(5,5,5,4) 


-25.55(05) 


-25.70(07) 


-25.80(06) 


-25.90(07) 


f 5 ,5 ,5 ,5 ) 


-25.29(05) 


-25.44(07) 


-25.54(06) 


-25.64(07) 


(6,5,5,5) 


-25.09(05) 


-25.25(06) 


-25.34(06) 


-25.45(07) 


(6,6,5,5) 


-24.90(04) 


-25.05(06) 


-25.15(06) 


-25.25(06) 


(6,6,6,5) 


-24.71(04) 


-24.86(06) 


-24.96(06) 


-25.06(06) 


(6,6,6,6) 


-24.52(04) 


-24.67(06) 


-24.77(06) 


-24.87(06) 


(7,6,6,6) 


-24.42(04) 


-24.57(06) 


-24.67(05) 


-24.77(06) 


(7,7,6,6) 


-24.31(04) 


-24.46(06) 


-24.56(05) 


-24.66(06) 


(7,7,7,6) 


-24.19(04) 


-24.34(06) 


-24.44(05) 


-24.54(06) 


(7,7,7,7) 


-24.07(04) 


-24.22(06) 


-24.32(05) 


-24.42(06) 



Table 5 

The vacuum polarisation Cn(y ,a?Ti) on the 16^ lattice at /3 = 6.0. 
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n 


n = 0.1345 


(1,0,0,0) 


-41.64(22) 


(1,1,0,0) 


-40.67(19) 


(1,1,1,0) 


-39.85(17) 


(1,1,1,1) 


-39.16(16) 


(2,1,1,1) 


-37.58(14) 


(2,2,1,1) 


-36.51(13) 


(2,2,2,1) 


-35.72(13) 


(2,2,2,2) 


-35.09(12) 


(3,2,2,2) 


-34.20(12) 


(3,3,2,2) 


-33.52(11) 


(3,3,3,2) 


-32.99(11) 




-32.54(11) 


(4,3,3,3) 


-31.97(11) 


(4,4,3,3) 


-31.50(10) 


(4,4,4,3) 


-31.11(10) 


(4,4,4,4) 


-30.77(10) 


(5,4,4,4) 


-30.36(10) 


(5,5,4,4) 


-30.01(10) 


(5,5,5,4) 


-29.70(10) 




-29.43( 9) 


(6,5,5,5) 


-29.11( 9) 


(6,6,5,5) 


-28.84( 9) 


(6,6,6,5) 


-28.59( 9) 


(6,6,6,6) 


-28.36( 9) 


(7,6,6,6) 


-28.11( 9) 


(7,7,6,6) 


-27.88( 9) 


(7,7,7,6) 


-27.68( 9) 


(7,7,7,7) 


-27.48( 9) 



Table 6 

The vacuum polarisation Cn(y ,a?Ti) on the 32^ lattice at (3 
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n 


n = 0.1344 


K = 0.1349 


K = 0.1352 


(1,0,0,0) 


-40.24(14) 


-40.87(17) 


-41.16(16) 


(1,1,0,0) 


-38.78(12) 


-39.17(14) 


-39.35(12) 


(1,1,1,0) 


-37.74(11) 


-38.03(13) 


-38.17(11) 


(1,1,1,1) 


-36.96(10) 


-37.19(12) 


-37.30(11) 


(2,1,1,1) 


-35.23(09) 


-35.40(10) 


-35.47(09) 


(2,2,1,1) 


-34.16(08) 


-34.31(09) 


-34.38(08) 


(2,2,2,1) 


-33.39(07) 


-33.54(08) 


-33.60(07) 


(2,2,2,2) 


-32.79(07) 


-32.94(08) 


-33.00(07) 


(3,2,2,2) 


-31.93(06) 


-32.08(07) 


-32.13(07) 


(3,3,2,2) 


-31.29(06) 


-31.44(07) 


-31.49(06) 


(3,3,3,2) 


-30.78(06) 


-30.94(07) 


-30.99(06) 




-30.35(06) 


-30.52(06) 


-30.57(06) 


(4,3,3,3) 


-29.81(06) 


-29.98(06) 


-30.03(06) 


(4,4,3,3) 


-29.37(06) 


-29.54(06) 


-29.59(06) 


(4,4,4,3) 


-28.99(05) 


-29.17(06) 


-29.22(06) 


(4,4,4,4) 


-28.67(05) 


-28.84(06) 


-28.90(05) 


(5,4,4,4) 


-28.29(05) 


-28.46(06) 


-28.52(05) 


(5,5,4,4) 


-27.96(05) 


-28.14(05) 


-28.20(05) 


(5,5,5,4) 


-27.67(05) 


-27.85(05) 


-27.91(05) 




-27.41(05) 


-27.58(05) 


-27.64(05) 


(6,5,5,5) 


-27.12(05) 


-27.30(05) 


-27.37(05) 


(6,6,5,5) 


-26.87(05) 


-27.05(05) 


-27.11(05) 


(6,6,6,5) 


-26.64(05) 


-26.82(05) 


-26.89(05) 


(6,6,6,6) 


-26.42(05) 


-26.60(05) 


-26.66(05) 


(7,6,6,6) 


-26.21(05) 


-26.39(05) 


-26.45(05) 


(7,7,6,6) 


-26.01(05) 


-26.19(05) 


-26.26(05) 


(7,7,7,6) 


-25.82(04) 


-26.00(05) 


-26.07(04) 


(7,7,7,7) 


-25.64(04) 


-25.82(05) 


-25.89(04) 



Table 7 

The vacuum polarisation Cn(y ,a?Ti) on the 24^ lattice at /? = 6.2. 
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n 


K = 0.1346 


n = 0.1350 


K = 0.1352 


(1,0,0,0) 


-41.88(16) 


-42.52(22) 


-42.84(18) 


(1,1,0,0) 


-40.38(15) 


-40.77(19) 


-40.97(15) 


(1,1,1,0) 


-39.33(14) 


-39.61(17) 


-39.75(14) 


(1,1,1,1) 


-38.53(13) 


-38.74(17) 


-38.86(14) 


(2,1,1,1) 


-36.82(11) 


-36.96(13) 


-37.02(11) 


(2,2,1,1) 


-35.73(10) 


-35.84(11) 


-35.89(10) 


(2,2,2,1) 


-34.93(09) 


-35.04(10) 


-35.09(09) 


(2,2,2,2) 


-34.31(08) 


-34.42(10) 


-34.46(08) 


(3,2,2,2) 


-33.44(08) 


-33.55(08) 


-33.58(08) 


(3,3,2,2) 


-32.79(07) 


-32.89(08) 


-32.93(07) 


(3,3,3,2) 


-32.27(07) 


-32.37(07) 


-32.41(07) 




-31.83(06) 


-31.94(07) 


-31.98(06) 


(4,3,3,3) 


-31.28(06) 


-31.39(06) 


-31.42(06) 


(4,4,3,3) 


-30.83(06) 


-30.94(06) 


-30.98(06) 


(4,4,4,3) 


-30.45(05) 


-30.56(06) 


-30.60(05) 


(4,4,4,4) 


-30.12(05) 


-30.23(06) 


-30.28(05) 


(5,4,4,4) 


-29.73(05) 


-29.84(05) 


-29.88(05) 


(5,5,4,4) 


-29.39(05) 


-29.51(05) 


-29.55(05) 


(5,5,5,4) 


-29.09(05) 


-29.21(05) 


-29.25(05) 




-28.83(04) 


-28.95(05) 


-28.99(05) 


(6,5,5,5) 


-28.53(04) 


-28.65(05) 


-28.69(04) 


(6,6,5,5) 


-28.26(04) 


-28.39(05) 


-28.43(04) 


(6,6,6,5) 


-28.02(04) 


-28.15(04) 


-28.19(04) 


(6,6,6,6) 


-27.80(04) 


-27.93(04) 


-27.97(04) 


(7,6,6,6) 


-27.56(04) 


-27.70(04) 


-27.74(04) 


(7,7,6,6) 


-27.35(04) 


-27.48(04) 


-27.52(04) 


(7,7,7,6) 


-27.15(04) 


-27.28(04) 


-27.32(04) 


(7,7,7,7) 


-26.96(04) 


-27.10(04) 


-27.14(04) 



Table 8 

The vacuum polarisation Cn(y ,a?Ti) on the 32^ lattice at /? = 6.4. 
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